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Abstract 

River  ice  velocity  measurements  are  fundamental  to  analyses  of  river  ice 
dynamics.  Ice  velocity  measurement  with  a  continuous-wave  Doppler  radar 
system  having  real-time  data  acquisition  and  digital  signal  processing  capa¬ 
bility  was  evaluated  during  a  river  breakup  and  a  frazil  run  on  the  Connecticut 
River.  This  system  can  be  rapidly  deployed,  requires  minimal  operator  interac¬ 
tion,  will  continuously  acquire,  process,  store,  and  display  ice  velocity  data  and 
does  not  depend  on  visibility  conditions.  In  parallel,  video  records  of  ice  motion 
were  obtained  atthe  same  location  for  later  manual  processing  and  comparison 
with  the  radar  results.  We  describe  the  Doppler  radar  system  and  obtain 
bounding  estimates  of  possible  measurement  errors.  The  principal  error  in 
Doppler  ice  velocity  measurement  is  due  to  the  beam  width  of  the  radar  antenna, 
and  an  analytical  method  is  developed  to  minimize  this  error.  Measured  ice 
velocities  ranged  from  1  to  2 .5  m/s  during  the  river  breakup  andfrom  0.5  to  0.65 
m/s  in  the  frazil  run.  Quantitative  comparisons  between  the  radar  and  video 
results  show  fundamental  agreement  between  these  measurement  methods, 
and  demonstrate  that  Doppler  radar  is  an  effective,  efficient,  and  precise  tool  for 
obtaining  river  ice  velocities  over  the  full  range  of  possible  ice  and  velocity 
conditions. 


Cover:  Frazil  ice  run  on  the  Connecticut  River,  28  December  1993.  Top:  From 
the  CornIsh-WIndsor  Bridge  radar  site  looking  upstream.  Bottom:  From 
the  teft  bank  video  site  looking  across  the  river.  (Photos  by  M.  Ferrick.) 


For  conversion  of  Si  units  to  non-SI  units  of  measurement  consult  Standard 
Practice  for  Use  of  the  International  System  of  Units  (SI),  ASTM  Standard  E380- 
93,  published  by  the  American  Societyfor  Testing  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 
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INTRODUCTION 

Accurate  ice  velocity  measurement  during  river 
breakup  in  spring  is  needed  for  analysis  of  river 
and  ice  dynamics  and  for  flood  hazard  or  water 
resource  assessment.  Ice  velocity  measurements 
presented  by  Ferrick  et  al  (1993)  demonstrated  the 
spatial  and  temporal  variability  of  the  motion  dur¬ 
ing  a  dynamic  breakup  of  the  Connecticut  River. 
Similar  though  generally  slower  and  more  grad¬ 
ually  variable  motion  occurs  during  the  freeze-up 
of  a  river.  Ferrick  et  al.  (1992)  described  the  acqui¬ 
sition  of  ice  velocity  data  by  locating  a  leiigth  refer¬ 
ence  on  the  ice,  recording  the  ice  motion  on  video¬ 
tape,  developing  a  computer-generated  reference 
grid  that  is  dubbed  over  the  videotape,  and  timing 
the  motion  of  individual  floes  through  the  grid. 
The  shortcomings  of  this  method  are  that  it  is 
laborious  and  it  requires  visual  observation  of  the 
ice,  which  can  be  precluded  by  fog  and  darkness. 
Ice  scour  of  the  bed  and  banks  frequently  prevents 
the  use  of  measurement  methods  that  require 
sensor  and  cable  placement  in  the  river.  Radar  sys¬ 
tems  have  been  applied  to  hydrological  investiga¬ 
tions,  including  the  measurement  and  mapping  of 
ice  thickness  on  rivers  and  lakes  (Arcone  and 
Delaney  1987,  Arcone  1991,  Yankielun  et  al.  1992, 
1993)  and  the  continuous  monitoring  of  river  stage 
(Yankielun  and  Ferrick  1993).  Doppler  radar  tech¬ 
niques  for  velocity  measurement  are  well  estab¬ 
lished  (Barton  1964,  Skolnik  1980)  and  provide  a 
remotely  mounted,  readily  installed,  continuously 
recording  alternative  for  ice  velocity  measurement 
that  is  not  affected  by  light,  fog,  or  ice  motion  con¬ 
ditions. 

A  continuous  wave  (CW)  Doppler  radar  with 
real-time  data  acquisition  and  digital  signal  pro¬ 
cessing  capability  was  moimted  on  the  Cornish 
(N.H.)-Wmdsor  (Vt.)  bridge  over  the  Connecticut 


River  during  the  March  1993  ice  breakup  and 
again  in  December  1993  during  a  frazil  ice  run 
prior  to  ice  cover  formation.  Ice  velocity  data 
were  continuously  acquired,  processed,  dis¬ 
played,  and  stored.  The  system  can  be  transport¬ 
ed  in  a  vehicle  and  rapidly  redeployed  to  permit 
data  acquisition  at  a  number  of  locations  during  a 
single  event  with  minimal  operator  interaction. 
Simultaneous  ice  velocity  data  were  obtained  at 
the  same  site  with  video  techniques  during  both 
of  these  events. 

In  this  paper  we  describe  the  Doppler  radar 
system,  analyze  the  contributions  to  measure¬ 
ment  error  inJierent  in  the  system,  determine  er¬ 
ror  bounds,  and  develop  data  processing  meth¬ 
ods  to  minimize  the  error  in  the  velocity  measure¬ 
ment.  Quantitative  comparisons  of  measured  ice 
velocity  are  made  between  the  radar  and  video 
techniques  for  both  the  breakup  and  frazil  events. 
Methods  for  enhancing  the  accuracy  of  Doppler 
radar  ice  velocity  measurement  and  its  ease  of 
use  follow  from  die  analysis. 

SYSTEM  DESCRIPTION 

The  primary  requirement  for  a  successful 
radar  measurement  is  that  the  target  provide  suf¬ 
ficient  backscatter  of  the  incoming  wave.  If  the 
surface  area  illuminated  by  the  radar  were  per¬ 
fectly  smooth,  the  incident  energy  would  be  re¬ 
flected  away  from  the  radar  antenna.  However, 
the  inherent  roughness  of  frazil,  sheet,  and  rub¬ 
ble  ice  provides  backscatter  towards  the  antenna 
(Fig.  1).  The  fundamental  relationship  between 
transmitted  and  received  power  of  a  radar  signal 
traveling  through  free  space  was  described  in 
Skolnik  (1980).  Lewis  et  al.  (1987)  discussed  the 
specific  problem  of  radar  detection  of  floating  ice 
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Figure  1.  The  effect  of  surface  roughness  on  radar  hackscatter,  and  defining 
terms  and  parameters. 


masses  on  large  bodies  of  water.  The  backscat- 
tered  power  received  Pj-  (W)  is  related  to  system 
and  target  parameters  as 

{AnfR^L' 

where  =  power  transmitted  (W) 

G  =  antenna  gain 
X  =  radar  wavelength  (m) 

A  =  area  illuminated  by  the  radar  (m^) 

R  =  range  to  the  target  (m) 

Go  =  normalized  cross  section  of  the  tar¬ 
get 

L  =  system  losses. 


An  ideal  reflector  returns  all  incident  energy  back 
to  the  source,  corresponding  to  a  Go  of  1.  The 
radar  cross  section  of  a  target  is  defined  as  the 
normalized  area  of  an  ideal  reflector  that  would 
return  the  same  signal  strength  as  the  target.  Sev¬ 
eral  factors  affect  the  radar  cross  section  of  a 
given  target,  including  transmitted  wavelength, 
target  geometry,  aspect,  and  reflectivity.  A  cali¬ 
brated  radar  system  is  necessary  to  measure  Gq 
for  various  target  conditions.  Moving  river  ice 
presents  an  inhomogeneous  target  with  complex 
and  variable  geometry,  and  the  radar  cross  sec¬ 
tion  can  vary  widely. 

An  additional  requirement  for  Doppler  veloc- 
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ity  measurement  is  that  the  frequency  shift 
obtained  must  be  resolvable  by  the  radar  system. 
In  a  CW  Doppler  system  a  single-frequency  RF 
carrier  is  transmitted  toward  the  target  at  an  obli¬ 
que  angle  (Fig.  1).  When  the  target  is  moving,  the 
backscattered  signal  exhibits  an  apparent  shift 
from  the  transmitted  carrier  frequency — ^higher  if 
movement  is  toward  the  radar  and  lower  if  it  is 
away  from  it.  A  low-level  sample  of  the  transmit¬ 
ted  carrier  frequency /s  is  used  as  a  reference  and 
mixed  with  the  backscattered  signal  received.  The 
mixing  process  results  in  four  frequencies:  the 
reference  signal,  the  backscattered  signal,  and  the 
sum  and  difference  of  these  frequencies.  The  dif¬ 
ference  frequency  resulting  from  the  mixing  is  the 
Doppler  frequency /d  proportional  to  the  velocity 
of  the  target.  The  velocity  v  (m/s)  of  a  moving  tar¬ 
get  is  obtained  as  (Roddy  1984) 

(2) 

2cosa  cosp 

where  X  =  wavelength  (m)  of  the  RF  source 

a  =  horizontal  off-stream  angle  as  shown 
in  Figure  1. 

P  =  vertical  look-angle  measured  down¬ 
ward  from  the  horizontal 

The  wavelength  and  frequency  of  the  source  are 
related  as  ^  =  c//s,  where  c  is  the  speed  of  light  in 
a  vacuum  (3  x  10^  m/s).  The  power  transmitted 
by  the  radar  is  not  localized  at  the  point  on  the 
target  surface  aligned  with  the  center  of  the  an¬ 
tenna  boresight.  Instead,  it  is  distributed  over  an 
area  on  the  surface,  with  the  incident  power  di¬ 
minishing  with  distance  from  the  boresight.  The 
3-dB  beam  width  shown  in  Figure  1  is  defined  as 
the  angular  width  of  the  radar  beam  at  the  half¬ 
power  points,  as  illustrated  in  Figure  2. 

The  Doppler  radar  implementation  used  for 
these  experiments  (Fig.  3)  includes  a  radar  front- 
end  and  a  data  acquisition  and  processing  com¬ 
puter.  The  radar  front-end  system  (expanded  in 
Fig.  3)  consists  of  a  standard  gain  horn  antenna 
that  is  connected  via  a  circulator  to  a  single-fre¬ 
quency,  continuous-wave  microwave  source  and 
a  single-ended  diode  mixer.  The  circulator  is  a  fer¬ 
romagnetic  device  that  directionally  commutates 
power  on  a  sequential  port-to-port  basis  and  per¬ 
mits  simultaneous  noninterfering  transmission 
and  reception  of  the  radar  signal  using  a  single 
antenna.  Ideally,  all  power  entering  port  1  is 
transmitted  to  port  2,  power  entering  port  2  is 
transmitted  to  port  3,  and  power  entering  port  3  is 
transmitted  to  port  1.  Port-to-port  isolation  in  the 


Figure  2.  The  3-dB  beam  width  of  a  radar  antenna. 


Data  Acquisition  and 
Processing  Hardware 


Figure  3.  Fundamental  components  of  the  Doppler 
radar  system  with  an  expanded  view  of  the  front-end. 


3 


reverse  direction  is  not  complete,  with  a  typical 
power  "leakage"  of  about  -20  dB.  The  leakage  be¬ 
tween  ports  1  and  3  is  used  as  the  reference  signal 
for  the  mixer.  The  front-end  assembly  is  mounted 
on  a  tripod  that  can  be  accurately  positioned  with 
an  inclinometer  to  within  ±1°  in  the  vertical  plane. 
Horizontally,  the  radar  is  visually  aimed  upstream 
into  the  flow,  and  alignment  is  adjusted  during  the 
ice  motion  by  manually  panning  the  antenna  until 
the  radar  response  indicates  a  velocity  maximum. 
Data  acquisition  and  display  were  performed 
with  a  33-MHz  80386  DOS  computer  system.  A 
DSP  card  acquires  16-bit  data  and  displays  the  dig¬ 
ital  signal  of  the  velocity  spectrum  in  real  time. 
Data  were  recorded  continuously  on  one  track  of  a 
four-track  digital  audio  tape  (DAT)  recorder  for 
later  playback,  processing,  and  analysis.  The  DAT 
recorder  was  also  used  to  record  concurrent  river 
stage  data  from  a  millimeter-wave  FM-CW  radar 
(Yankielun  and  Ferrick  1993),  event  timing,  and  a 
voice  channel  for  a  descriptive  narrative. 

After  completion  of  a  survey,  the  raw  Doppler 
data  were  processed  and  displayed.  Each  radar 
scan  was  digitized  to  provide  1024  time  series 
samples,  transformed  into  a  power  spectrum  and 
processed  with  a  Hanning  window  to  suppress 
the  effect  of  spectral  sidelobes  that  could  mask 
lower-level  signals.  The  processed  power  spec¬ 
trum  can  be  displayed  in  either  a  single-scan  for¬ 
mat  or  as  a  continuous  series  of  scans  in  spectro- 
graphic  form  (Fig.  4).  In  a  spectrogram,  discrete 


Figure  4.  Power  spectra  of  single  and  multiple  scans 
that  indicate  the  relationship  between  the  amplitude 
threshold  setting  and  the  bandwidth  of  Doppler  fre¬ 
quencies. 


signal  magnitude  quanta  are  represented  by  a 
range  of  color  or  gray  scale.  With  16-shade  gray 
scale  graphics,  the  maximum  signal  magnitudes 
appear  black,  and  intermediate  levels  appear  as 
lighter  shades  of  gray.  Below  a  preset  magnitude 
threshold  all  is  shown  as  white.  The  levels  can  be 
set  in  the  DSP  software  to  display  the  Doppler  fre¬ 
quency  clearly.  A  multicolor  spectrographic  dis¬ 
play  provides  a  256-shade  color  gradient  that  sig¬ 
nificantly  improves  graphical  resolution. 

ERROR  ANALYSIS 

To  assess  the  capability  of  Doppler  radar  to 
measure  ice  velocity  we  must  identify  and  quanti¬ 
fy  the  sources  of  error  inherent  in  the  method  and 
minimize  these  errors  for  minimum  total  error.  It 
would  be  useful  to  combine  the  individual  errors 
and  obtain  an  upper  bound  on  the  total  error.  We 
take  i;  in  eq  2  as  the  dependent  variable,  and  write 
the  total  differential  dv  as 

An  upper  bound  on  the  total  error  is  obtained  as 
the  sum  of  the  absolute  values  of  the  terms  in  eq  3, 
when  each  term  represents  an  individual  upper 
bound.  The  differentials  of  the  independent  vari¬ 
ables  will  be  replaced  by  finite  quantities  that  we 
assume  are  small  enough  for  eq  3  to  provide  an 
accurate  estimate  of  each  component  of  error* 

The  vertical  angle  and  horizontal  angle  must 
be  known  precisely  to  obtain  an  accurate  velocity 
measurement.  The  partial  derivatives  of  velocity 
in  eq  2  with  respect  to  a  and  (5  have  the  same  form: 

•^  =  i;tanp;  -^  =  i?tana.  (4) 

9p  da 

When  maximum  error  is  being  evaluated,  the  val¬ 
ues  of  a  and  P  in  eq  4  should  be  the  most  probable 
angle  plus  the  estimate  of  angular  error.  The 
wavelength  of  the  CW  RF  source  is  known  and 
should  be  precise  and  stable  for  accurate  mea¬ 
surement.  Wavelength  error  is  the  difference  be¬ 
tween  the  actual  and  the  measured  radar  source 
carrier  wavelength.  The  quantized  nature  of  digi¬ 
tal  sampling  is  a  source  of  error  in  the  Doppler  fre¬ 
quency.  The  partial  derivative  of  velocity  in  eq  2 
with  respect  to  wavelength  X  and  Doppler  fre¬ 
quency/d  have  the  same  form: 
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We  can  now  rewrite  eq  3  as  an  upper-bound 
percentage  error  in  velocity,  expressing  each  of 
the  differentials  as  finite  errors  A: 

—  =  Itan  p  •  Ap|  +  |tana  •  Aa|  +  —  +  .  (6) 

Our  vertical  angle  measurement  by  inclinometer 
of  25®  was  accurate  to  within  ±2®,  and  the  hori¬ 
zontal  angle  alignment  to  the  maximum  velocity 
at  0®  was  accurate  to  about  ±5®.  The  upper-boimd 
percentage  error  contributions  of  these  terms  are 
then  0.018  and  0.0076,  respectively.  A  source  fre¬ 
quency  of  5  GHz  and  an  upper  bound  frequency 
error  of  100  MHz  yield  a  percentage  error  in  X,  of 
0.02.  The  percentage  error  in/d  is 


^=-1  (7) 

/d  ^ 

where  N  is  the  ntimber  of  samples  acquired  and 
processed  per  second  for  each  Doppler  frequency 
determination.  In  our  system,  the  sampling  fre¬ 
quency  N  was  1024  Hz,  corresponding  to  a  negli¬ 
gible  digital  sampling  error  of  less  than  0.001. 
The  sum  of  the  absolute  values  of  the  upper- 
boimd  error  estimates  in  eq  6  yields  I  Au/ul  equal 
to  0.047. 

The  beam  width  of  the  radar  antenna  intro¬ 
duces  another  potential  error,  in  addition  to  those 
given  by  eq  6.  The  antenna  beam  width  corre¬ 
sponds  to  an  area  on  the  surface  of  the  ice  that  is 
illuminated  by  the  radar  and  returns  backscatter 
to  the  antenna.  Because is  angle-dependent,  a 
band  of  Doppler  frequencies  represents  the  back¬ 
scatter  from  across  the  illuminated  area.  The  near 
edge  of  the  beam  corresponds  to  the  largest  verti¬ 
cal  angle  and  the  upper-boimd  error.  Following 
the  approach  taken  above,  we  obtain  an  upper- 
boimd  estimate  of  this  error  by  replacing  p  by  p  + 
0/2  in  eq  2  and  taking  the  partial  derivative  of 
velocity  with  respect  to  0/2: 


0"| 

2  j  3(6/ 2) 


=  »tan 


(8) 


When  is  6  greater  than  10°,  the  beam-width  error 
given  by  eq  8  is  potentially  larger  than  the  sum  of 
the  upper-bound  errors  in  eq  6. 


VELOCITY  FROM  A  DOPPLER  FREQUENCY 
BAND  AND  A  VIDEO  RECORD 

Systems  with  large  beam  width  antermas  re¬ 
quire  additional  data  analysis  to  obtain  an  accu¬ 


rate  velocity  from  the  Doppler  frequency  band  at 
each  selected  time  throughout  the  record.  If  the 
actual  vertical  angle  P  representing  any  point  in 
the  frequency  band  were  known,  the  velocity  at 
the  corresponding  time  could  be  obtained  from  eq 
2.  However,  these  angles  are  not  generally  known. 
The  apparent  far  and  near  edge  returns  are  the 
most  readily  identifiable  locations  on  the  frequen¬ 
cy  band.  The  width  of  this  band  depends  on  the 
antenna  beam  width  and  the  threshold  setting  on 
the  power  spectrum  of  the  backscattered  signal. 
Our  initial  step  in  data  reduction  is  to  replace  the 
frequency  scale  of  the  data  band  by  a  velocity 
scale  using  eq  2  with  the  boresight  angle  as  p.  The 
center  of  ttiis  band  is  the  approximate  ice  velocity, 
and  it  can  be  obtained  immediately.  For  more  pre¬ 
cise  work  during  postprocessing,  the  apparent  far 
and  near  edge  velocities  and  v_,  corresponding 
to  the  upper  and  lower  edges  of  the  band,  are  re¬ 
lated  to  each  other  and  the  desired  velocity  v  as 

COSP  COSp  (Q\ 

v= — — 7 — 
cosj^P-^j  cos[p  +  ^J 

where  6a  is  an  apparent  beam-width  angle.  The 
velocity  v  corresponding  to  the  anteima  boresight 
is  bounded  above  by  and  below  by  v_.  We  ob¬ 
tain  the  difference  between  the  edge  velocities  as 
a  percentage  of  the  velocity  v  from  eq  9  as 


V  COsP 


(10) 


The  near  and  far  edge  velocity  ratios  vj v  and  v+/ 
V  are  plotted  in  Figure  5  as  a  function  of  0a  for  a 
range  of  boresight  angles  p.  The  potential  beam- 


ea(°) 

Figure  5.  Relationship  between  the  velocity  ratio  and 
%for  both  the  near  and  far  edges  of  the  beam  width 
and  a  range  of  p  values. 
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Figure  6,  Example  that  quantifies  the  near-  and  far-edge  velocity  ratio  with  the  boresight 
velocity  for  a  specific  case,  p  =  25°  and  0^  =  20  °,  and  relates  them  to  specific  locations  on  the 
spectrogram. 


width  errors  increase  significantly  with  for  all 
boresight  angles  and  increase  with  p  for  all  beam 
widths.  A  sample  Doppler  velocity  spectrogram 
is  related  to  the  velocity  ratio  relationship  for  p  = 
25°  and  0a  =  20°  in  Figure  6. 

Manipulating  the  relationship  between  the 
edge  velocities  in  eq  9  we  solve  for  0a  as 


0 


a 


=  2tan  ^ 


-t  +  1 


r- 


J 


tanp 


(11) 


If  both  upper  and  lower  bounding  traces  are 
drawn  consistently,  0a  should  vary  randomly 
within  a  narrow  range.  Velocity  estimates  Vi  and 
V2  are  obtained  with  the  mean  0a  as 


cosp 


cos 


^  2 


Vj  =- 


cosp 


cos|  p  +  -^ 


(12) 


Finally,  the  measured  Doppler  velocity  v,  repre¬ 
senting  the  bulk  ice  velocity  within  the  radar 
footprint,  is  taken  as  the  average  of  Vi  and  i?2‘ 

A  reference  grid  for  video  ice  velocity  mea¬ 
surement  was  established  by  placing  visual  tar¬ 
gets  on  the  winter  ice  cover  at  known  spacing. 
The  video  records  of  the  breakup  and  frazil  run 
events  were  obtained  by  timing  ice  floes  through 
this  grid.  As  in  previous  work,  we  assume  that 
the  measured  data  contain  a  slowly  varying  com¬ 


ponent  that  provides  the  global  variation  of  ice 
velocity  with  time  and  a  comparatively  rapid 
variation  of  small  amplitude  that  represents 
measurement  error  or  noise  in  the  data.  Separat¬ 
ing  these  components  requires  the  development 
of  a  function  that  contains  most  of  the  informa¬ 
tion  in  the  data  and  little  of  the  error  or  noise. 
The  video  data  were  fitted  with  orthogonal  poly¬ 
nomials  of  increasing  order  until  the  maximum 
differences  ceased  to  diminish  and  the  time  se¬ 
ries  of  these  differences  appeared  random.  The 
polynomial,  as  a  single-valued  function  of  time, 
permits  quantitative  comparison  of  the  radar 
and  video  records. 

Prowse  and  Demuth  (1991)  emphasized  the 
importance  of  a  viewing  position  with  large  free¬ 
board  and  accurate  vertical  angle  measurement 
to  obtain  accurate  ice  velocity  measurement  with 
a  related  optical  method.  The  water  levels  during 
breakup  were  much  higher  than  those  during 
target  setup,  occupying  30%  of  the  setup  free¬ 
board.  The  resulting  change  in  vertical  angle 
from  the  camera  to  all  points  on  the  ice  caused 
distortion  of  the  original  grid  by  length-scale 
changes  that  depend  on  position  in  the  cross  sec¬ 
tion  of  the  river.  The  open  water  during  the  frazil 
ice  nm  was  at  a  level  comparable  to  that  of  the 
winter  ice  cover.  The  data  from  both  events  were 
obtained  along  a  single  streamline  in  the  river, 
providing  consistent  results  that  can  be  corrected 
for  a  change  in  length  scale.  The  velocity  data  ob¬ 
tained  from  the  original  grid  v^^id  are  related  to 
the  actual  velocity  v  as 


6 


^Xgrid  ^X  +  AX 
^8^**  Af  At 


=  v  +  Av, 


(13) 


where  X  is  the  actual  length  scale.  At  is  the  transit 
time  of  a  floe  over  X,  and  Av  is  a  velocity  error 
corresponding  to  length-scale  error  AX.  The 
actual  video  velocity  v  can  then  be  obtained  as 


^grid  ^^grid  ^ 


(14) 


where  C  =  1  indicates  no  grid  distortion,  and 
increasing  grid  error  is  indicated  as  C  departs 
from  1. 


Ice  breakup — March  1993 

A  30-minute  Doppler  radar  ice  velocity  record 
and  a  simultaneous  video  ice  velocity  record 
were  obtained  for  the  Connecticut  River  at  the 
Cornish-Windsor  bridge  on  30  March  1993. 
Specifications  that  characterize  the  radar  system 
used  are  given  in  Table  1.  The  radar  antenna  was 
moxmted  about  6  m  above  the  ice  and  aimed  di¬ 
rectly  upstream  with  an  angle  p  of  25°  below  hor¬ 
izontal.  Concurrently,  a  video  camera  was 
mounted  at  the  top  of  the  left  bank  and  aimed 
normal  to  the  direction  of  the  flow.  The  stream¬ 
line  where  video  data  were  obtained  passed 
through  the  footprint  of  the  radar.  Ice  breakup 


Table  1.  Doppler  system  configurations. 

_ Event _ 


Breakup  Frazil  run 


Radar  band 

C 

Ka 

Frequency 

5  GHz 

30  GHz 

Wavelength 

6  cm 

1  cm 

Antenna  gain 

15  dB 

27  dB 

3  dB  antenna  beam  width 

20° 

9° 

Microwave  power 

20  dBm 

20  dBm 

was  underway  at  a  river  flow  of  approximately 
850  mVs  when  measurements  began.  For  most 
of  this  record,  the  ice  was  densely  packed  rubble 
that  provided  strong  radar  reflections  and  many 
video  targets.  Toward  the  end  of  the  sampling 
period,  the  surface  ice  concentration  diminished 
and  the  ice  velocity  increased. 

The  ice  velocity-time  trace  obtained  from  the 
Doppler  frequency  band  is  presented  in  Figure  7. 
Digital  signal  processing  software  identified  the 
peak  return  from  each  radar  scan  in  this  record. 
These  peak  detector  data,  representing  the  high¬ 
est  magnitude  power  reflections  from  within  the 
footprint  illuminated  by  the  radar,  are  also  plot¬ 
ted  in  Figure  7.  Both  of  these  data  sets  were 
processed  by  first  identifying  the  upper  v+  and 
lower  v_  edges  of  the  trace  at  20-  or  30-s  intervals 
throughout  the  record;  the  results  are  shown  in 
Figure  7.  The  peak  returns  generally  originated 
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Figure  7.  Data  band  of  raw  Doppler  velocities  and  peak  returns  for  the  ice  breakup  event  and 
boundaries  that  were  drawn  for  each  of  these  data  sets. 
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Figure  8.  Calculated  ©«  as  a  function  of  time 
for  the  raw  and  peak  radar  data  records  for  the 
breakup  event.  These  values  result  from  an 
equal  boresight  velocity  requirement  for  the 
near  and  far  edge  traces. 

close  to  the  near  edge  of  the  footprint,  consistent 
with  the  inverse  relationship  between  backscat- 
tered  power  and  range  to  the  fourth  power  given 
in  eq  1. 

Using  eq  11  we  calculated  Og  for  both  the  raw 
and  peak  detector  data  records;  we  present  the 
respective  results  in  Figure  8  as  a  function  of 
time.  For  the  raw  data  record  the  mean,  median, 
and  mode  of  %  are  approximately  equal,  with  a 
mean  of  44.6°,  a  standard  deviation  of  2.2°,  and  a 
range  of  10.0°.  The  mean,  median,  and  mode  of 
9a  for  the  peak  data  record  are  also  approximate¬ 
ly  equal,  with  a  mean  of  26.0°,  a  standard  devia¬ 
tion  of  2.1°,  and  a  range  of  9.4°.  The  peak-finding 
algorithm  significantly  reduced  the  apparent 
beam  width  of  the  radar  without  altering  the  0a 


distribution.  The  V\  and  U2  results  and  the  differ¬ 
ence  between  these  values  nondimensionalized 
by  the  corresponding  mean  velocity  are  given  in 
Figure  9.  Except  for  a  few  points,  the  absolute 
value  of  this  dimensionless  difference  was  less 
than  0.04  in  both  records,  corresponding  to  an 
absolute  value  of  the  difference  between  v  and 
either  Vi  or  V2  of  less  than  0.02  v.  The  delineation 
of  the  limits  of  each  radar  trace  and  calculation  of 
a  mean  0  has  reduced  the  large  potential  beam- 
width  error  to  the  same  order  as  the  other  minor 
errors  inherent  in  Doppler  velocity  measure¬ 
ment.  The  peak  data  are  smooth  relative  to  the 
raw  data,  but  they  are  shifted  systematically  to 
lower  velocities,  following  the  downward  shift 
of  the  far  edge  of  the  trace.  The  resulting  mean 
velocity  of  the  peak  radar  record  is  1.35  m/ s,  and 
that  of  the  raw  radar  record  is  1.61  m/ s. 

The  velocity  data  from  the  video  record  are 
shown  with  the  selected  polynomial  fit  and  the 
corresponding  normalized  error  in  Figure  10. 
The  mean  and  median  of  the  normalized  error 
were  very  close  to  zero,  with  a  standard  devia¬ 
tion  of  0.026.  The  maximum  normalized  error  of 
about  0.06  V  is  within  the  measurement  error 
boimd  for  this  method  (Ferrick  et  al.  1992).  Com¬ 
parison  of  the  radar  and  video  velocities  reveals 
that  the  video  results  are  consistently  low.  For 
ease  of  comparison  we  assume  that  this  system¬ 
atic  difference  is  due  entirely  to  video  grid  dis¬ 
tortion.  We  obtain  C  in  eq  14  from  the  ratio  of  the 
raw  radar  to  video  polynomial  mean  velocity 
values  as  1.067. 


Figure  9.  Ice  breakup  boresight  velocities  and  their  normalized  difference  obtained  from  the 
edge  traces,  xvilh  the  corresponding  mean  0^  of  both  the  raw  and  peak  radar  data  records. 
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Figure  10.  Ice  velocity  data,  best  polynomial  fit, 
and  normalized  difference  between  the  data  and 
polynomial  of  the  video  ice  breakup  record. 


Figure  11.  Comparison  of  ice  velocity  results 
for  breakup  after  multiplying  the  video  and 
peak  radar  results  by  a  constant  that  equates 
the  mean  velocity  of  each  record. 


Almost  identical  correlation  coefficients  were 
obtained  for  all  pairs  of  the  raw  radar,  peak  radar, 
and  polynomial  video  velocity  records,  ranging 
between  0.933  and  0.950.  The  ice  velocity-time 
traces  obtained  with  each  method  are  superposed 
in  Figure  11,  following  multiplication  of  the  peak 
return  results  by  1.191  and  the  video  polynomial 
results  by  1.067  to  correct  for  the  mean  offset  from 
the  raw  radar  results.  The  root-mean-square 
(RMS)  difference  obtained  between  any  pair  of 
records  was  less  than  7.7  cm/s,  corresponding  to 
dimensionless  RMS  velocity  differences  of  0.033, 
0.037,  and  0.048  between  the  video-raw,  video¬ 


peak,  and  raw-peak  records,  respectively  Given 
the  various  measurement  errors,  these  compari¬ 
sons  suggest  fundamental  agreement  of  the 
methods.  However,  additional  development  of 
the  peak  return  data  reduction  method  is  needed 
because  of  the  significant  imderprediction  of  the 
mean  velocity. 

Ice  velocities  exceeding  5  m/s  during  breakup 
have  been  reported  on  some  rivers.  The  maxi¬ 
mum  Doppler  frequency  that  can  be  obtained 
with  a  given  system  is  limited  by  the  Nyquist 
sampling  rate,/^  <  N/2.  With  our  system,  veloci¬ 
ties  of  up  to  6  m/s  can  be  obtained  with  Doppler 
frequencies  less  than  200  Hz,  indicating  that  the 
Nyquist  rate  does  not  restrict  velocity  measure¬ 
ment  at  the  high  end  of  the  observed  range.  Un¬ 
like  optical  methods,  the  large  changes  in  river 
stage  that  can  occur  during  breakup  do  not  affect 
Doppler  radar  velocity  measurement. 

Frazil  ice  run — December  1993 

Doppler  radar  ice  velocity  and  simultaneous 
video  ice  velocity  records  were  obtained  for  the 
Connecticut  River  at  the  Cornish-Windsor 
bridge  on  28  December  1993.  Frazil  pans  and 
floes  were  moving  downstream  during  the 
measurement  period  at  a  river  flow  of  about  170 
m^/s,  just  prior  to  ice  cover  formation.  Over  sev¬ 
eral  hours  of  data  collection,  the  river  flow  and 
ice  velocity  conditions  were  relatively  steady, 
and  we  selected  a  typical  1000-s  record  for  analy¬ 
sis.  The  video  setup  and  grid  were  the  same  as  in 
the  previous  event.  The  radar  antenna  was 
mounted  at  the  same  height  and  location  as  be¬ 
fore,  about  11  m  above  the  lower  water  surface. 
To  reduce  the  data  processing  requirements  of 
the  method,  an  antenna  with  a  much  narrower 
beam  width  was  used.  With  slower  ice  motion 
expected,  the  wavelength  of  the  radar  was  de¬ 
creased  to  improve  the  velocity  resolution.  Using 
eq  2,  the  Doppler  frequency  shift  needed  to  re¬ 
solve  velocity  differences  of  0.05  m/ s  is  only  1.5 
Hz  for  a  5-GHz  microwave  source.  The  corre¬ 
sponding  frequency  shift  increases  to  9  Hz  with  a 
30-GHz  millkneter-wave  source.  A  millimeter- 
wave  source  provides  sufficient  resolution  for  ice 
velocity  measurement  at  speeds  above  about  0.1 
m/s.  The  modified  radar  system  specifications 
are  given  in  Table  1. 

The  video  velocity  data,  best  fifth-order  poly¬ 
nomial  fit,  and  difference  between  the  data  and 
the  polynomial  normalized  by  the  polynomial 
velocity  are  given  in  Figure  12.  The  mean  video 
velocity  was  0.518  m/s,  and  the  normalized  dif- 
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Figure  12.  Video  and  radar  data,  best  polynomial  fits,  and  normalized  difference  between  the 
data  and  corresponding  polynomial  for  the  frazil  ice  run. 


ference  between  the  data  and  the  polynomial  was 
generally  less  than  0.06.  The  mean  and  median 
difference  are  both  zero,  with  an  RMS  difference 
of  0.022  m/s  and  a  maximum  difference  of  0.065 
m/s.  The  amplitudes  of  the  velocity  differences 
are  within  the  boimds  of  experimental  error,  and 
their  structure  is  random.  We  conclude  that  the 
polynomial  adequately  describes  the  video  data. 

The  full  width  of  the  Doppler  frequency  band 
was  a  difference  in  velocity  of  0.05  m/ s,  or  a  nor¬ 
malized  velocity  difference  of  0.09.  Data  repre¬ 
senting  the  midpoint  of  the  band  obtained  at  10-s 
intervals,  the  best  fourth-order  polynomial  fit  to 
these  data,  and  the  normalized  difference  be¬ 
tween  the  data  and  the  polynomial  are  presented 
in  Figure  12,  The  mean  radar  velocity  was  0.539 
m/s,  and  the  normalized  difference  was  general¬ 
ly  less  than  0.07,  The  mean  and  median  differ¬ 
ences  are  both  zero,  with  an  RMS  difference  of 
0.023  m/s  and  a  maximum  difference  of  0.067  m/ s. 
The  amplitudes  of  the  velocity  differences  are 
somewhat  larger  than  would  be  expected  from 
experimental  error,  and  their  structure  appears 
periodic.  However,  the  oscillations  in  this  record 
are  not  supported  by  the  video  data.  The  oscillat¬ 
ing  radar  velocities  occurred  when  the  open  wa¬ 
ter  area  in  the  radar  footprint  increased  and  the 
number  of  targets  decreased.  When  a  single  floe 
traverses  the  footprint  there  is  an  apparent 
change  in  velocity  from  high  to  low  caused  by  the 
relationship  between  vertical  angle  and  Doppler 
frequency.  Figure  13  shows  part  of  the  breakup 


Figure  13.  The  part  of  the  radar  velocity  record 
for  breakup  with  dispersed  icefloes  in  the  river. 
Note  the  movement  in  time  of  strong  returns  from 
higher  to  lower  velocity,  corresponding  to  the  move¬ 
ment  offices  through  the  antenna  footprint. 

record  during  the  period  of  dispersed  ice  motion. 
Bands  of  velocity  through  time  from  floes  produc¬ 
ing  strong  backscatter  are  clearly  visible. 

Small  underestimates  of  a  and  P  would  cause  a 
systematically  high  radar  velocity,  corresponding 
to  most  of  the  difference  between  the  Doppler 
and  video  mean  velocities.  In  addition,  strong  in¬ 
dividual  reflectors  that  passed  through  the  radar 
footprint  at  constant  speed  spent  almost  60%  of 
this  time  on  the  high-velocity  side  of  the  bore- 
sight,  potentially  introducing  a  velocity  bias. 
However,  for  consistency  and  ease  of  comparison 
we  again  assume  that  the  systematic  difference 
between  the  radar  and  video  velocities  is  due  to 
video  grid  distortion.  The  video  and  radar  veloc- 
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Figure  14.  Comparison  of  ice  velocity  results  for  the 
frazil  run  after  the  video  results  are  multiplied  by  a  con¬ 
stant  that  equates  the  mean  velocity  of  each  record. 

ity  polynomials  are  compared  in  Figure  14  after 
the  video  velocity  polynomial  is  multiplied  by 
1.041  to  correct  for  the  difference  in  the  means. 
Discounting  the  tail  of  the  radar  polynomial  that 
is  an  end-effect  of  the  final  oscillation,  the  RMS 
difference  between  these  curves  is  0.014  m/ s,  or 
0.026  when  normalized  by  the  mean  velocity.  Ex¬ 
cept  for  the  oscillations  induced  by  strong  and 
sparse  reflectors,  the  Doppler  radar  and  video 
data  agreement  is  well  wi^in  the  measurement 
errors  of  the  methods.  The  Doppler  radar  system 
again  provided  an  excellent  velocity  measure¬ 
ment  with  minimal  processing. 

CONCLUSIONS 

Doppler  ice  velocity  measurements  are  made 
remotely,  can  be  displayed  in  real  time,  and  are 
not  affected  by  light,  fog,  or  ice  motion  conditions. 
Doppler  radar  can  measure  and  resolve  the  veloc¬ 
ity  of  ice  moving  in  a  river  with  precision  com¬ 
parable  to  or  better  than  analysis-intensive  video 
techniques,  over  the  complete  range  of  ice  and  vel¬ 
ocity  conditions.  The  maximum  error  in  Doppler 
radar  measurement  of  river  ice  motion  is  about  5% 
of  the  velocity.  Agreement  between  the  mean 
Doppler  and  video  velocities  was  obtained  within 
7%  for  the  ice  breakup  case  study  and  4%  for  the 
frazil  run.  Video  grid  distortion  during  breakup 
and  sparse  radar  targets  during  the  frazil  run  are 
likely  causes  of  these  systematic  differences  in  vel¬ 
ocity.  The  RMS  differences  between  the  mean- 
corrected  velocity  curves  were  about  3%  for  both 
cases.  A  high  gain,  narrow-beam  antenna  im¬ 
proved  signal-to-noise  performance  of  our  Dop¬ 
pler  radar  system  and  minimized  data  processing 
requirements.  An  increase  in  the  source  frequency 
of  the  system  provided  a  proportional  increase  in 
the  velocity  resolution. 

Significant  cost  reduction  and  miniaturization 


can  now  be  realized  through  application  of  off- 
the-shelf  10-GHz  and  24-GHz  microwave  trans¬ 
ceiver  modules  consisting  of  a  horn  antenna,  cir¬ 
culator,  diode  mixer,  and  Gunn  oscillator.  These 
modules  only  require  connections  to  dc  power 
and  to  a  data  acquisition  computer  to  realize  a 
Doppler  radar  system,  encouraging  long-term  or 
permanent  mounting  of  sensors  at  strategic  points 
on  a  river  to  provide  early  warning  and  continu¬ 
ous  monitoring  of  ice  freeze-up  or  breakup.  A  lap¬ 
top  computer  equipped  with  a  PCMCIA  data  ac¬ 
quisition  card  enhances  system  portability  and  re¬ 
duces  power  requirements. 
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